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Abstract: Normal neutrosophic numbers (NNNs) are a significant tool of describing the 
incompleteness, indeterminacy and inconsistency of the decision-making information. In this paper, we 
firstly propose the definition and the properties of the NNNs, and the accuracy function, the score 
function and the operational laws of the NNNs are developed. Then, some operators are presented, 
including the normal neutrosophic Bonferroni mean (NNBM) operator, the normal neutrosophic 
weighted Bonferroni mean (NNWBM) operator, the normal neutrosophic geometric Bonferroni mean 
(NNGBM) operator, the normal neutrosophic weighted geometric Bonferroni mean (NNWGBM) 
operator. We also study their properties and special cases. Further, we put forward a multiple attribute 
decision making method which is based on the NNWBM and NNWGBM operators. Finally, an 
illustrative example is given to verify the practicality and validity of the proposed method. 

Keywords: Multiple attribute decision making; Normal neutrosophic numbers; Normal neutrosophic 
Bonferroni mean aggregation operator. 

1. Introduction 

As an important research branch of decision theory. Multiple attribute decision making (MADM) 
has a wide application in many areas. The multiple attribute decision making was firstly proposed and 
applied to select the investment policy of the enterprises by Churchman et al. [1]. However, because the 
fuzziness and indeterminacy of the information in real decision making are a common phenomenon, 
numerical values are inadequate or insufficient to model real-life decision problems. In same 
occasions, it can be more accurate to use the fuzzy numbers to describe the attribute values in fuzzy 
environment. Because of the fuzzy, indeterminate information in most of the multiple attribute decision 
making problems, Zadeh [2] firstly proposed the fuzzy set (FS) and used it to describe the fuzzy 
information. Atanassov [3] further presented the intuitionistic fuzzy set (IFS) where he added the 
non-membership function to the FS. Now, IFS has attracted more and more attentions in MADM and 
multiple attribute group decision making (MAGDM). However, due to the fact that the membership 
and non-membership function in IFS can only be shown by the accurate numbers, IFS is not 
appropriate to be used to solve the F-MADM problems. Atanassov [4], Gargov and Atanassov [5] 
further proposed the interval-valued intuitionistic fuzzy set (IVIFS). Liu and Zhang [6] defined 
triangular intuitionistic fuzzy numbers (TIFNs). Wang [7] proposed the trapezoidal intuitionistic fuzzy 
numbers (TrIFNs) and the interval trapezoidal intuitionistic fuzzy numbers (ITrIFNs). But, these 
extensions of IFS can just consider the indeterminacy. 

Smarandache [8] proposed the neutrosophic set (NS) in which an independent 
indeterminacy-membership function was added. In NS, the truth-membership function is the same 
function as the membership of IFS; likewise, the falsity-membership is the same function as the 
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non-membership of IFS. The indeterminacy-membership function is the pivotal difference between NS 
and IFS. The three parts, including the truth-membership, the indeterminacy membership, and the 
false-membership are completely independent. Based on this, Wang et al. [9] expressed the values of 
indeterminacy-membership, truth-membership and false-membership with the interval numbers, which 
has been defined the interval neutrosophic set (INS). On the basis of the Hamming and Euclidean 
distances. Ye [10] defined some similarity measures between INS. Then, a method about the multiple 
attribute decision making based on the similarity degree was proposed. 

In real-life world, the normal distribution is widely applied to a lot of fields. But both the IFS and 
INS cannot consider the normal distribution, so the researches about the normal fuzzy information are 
attracting more and more attentions. Yang and Ko[ll] firstly defined the normal fuzzy numbers(NFNs) 
to express the normal distribution phenomena. NFNs are more reasonable and realistic to express the 
decision-making information than FS in a random environment. Based on the NFNs and IFS, Wang and 
Li [12] proposed the normal intuitionistic fuzzy numbers (NIFNs) and defined its corresponding 
operations, the stability factor, the score function and so on. However, there have not been researches 
about the combination NFNs with NNs. 

Now, more and more researchers pay attention to the information aggregation operators, which 
have become a popular research topic [26-30]. Bonferroni [13] firstly proposed the Bonferroni mean 
(BM) operator which can catch the interrelationship between the input arguments, BM has been applied 
in many application domains and attracted more and more attentions from the researchers. Yager [14] 
proposed some generalizations about the BM, such as the ordered weighted averaging (OWA) operator 
[15] and Choquet integral [16]. Yager [17] and Beliakov et al. [18] defined another generalized form of 
BM. Nevertheless, Zhu et al. [19] proposed the geometric Bonferroni mean (GBM) in which both the 
BM and geometric mean (GM) are considered. 

Up to now, there is no research on the normal neutrosophic decision-making problems considering 
the interrelationship between the input normal neutrosophic arguments. Therefore, it is necessary to 
pay more attention to this issue. Because the BM operator can consider the interrelationship between 
the attributes, and the NNNs have the advantages of considering the normal random information and 
the neutrosophic variables, which can handle the incomplete, inconsistent and indeterminate 
information. In this paper, we extend the Bonferroni mean to aggregate the normal neutrosophic 
variables by combining BM aggregation operator with NNNs. We firstly propose two aggregation 
operators called the normal neutrosophic Bonferroni mean (NNBM) operator and the normal 
neutrosophic geometric Bonferroni mean (NNGBM) operator for aggregating the normal neutrosophic 
numbers. Then we study some properties of them and discuss their some special cases. For the 
situations in which the input arguments have different weight, we then develop the normal 
neutrosophic weighted Bonferroni mean (NNWBM) operator and the normal neutrosophic weighted 
geometric Bonferroni mean (NNWGBM) operator, and then, we propose two procedures for multiple 
attribute decision making under the environments of the NNNs based on the proposed operators. 

The remainder of this paper is constructed as follows. In the next section, we introduce some basic 
concepts of the NNNs, some operational laws and the prominent characteristics of NNNs. In section 3, 
some aggregation operators on the basis of the normal neutrosophic numbers are proposed, such as the 
normal neutrosophic Bonferroni mean (NNBM) operator, the normal neutrosophic weighted 
Bonferroni mean (NNWBM) operator, the normal neutrosophic geometric Bonferroni mean (NNGBM) 
operator, the normal neutrosophic weighted geometric Bonferroni mean (NNWGBM) operator, and 
their properties are discussed. In section 4, a multiple attribute decision making method on the basis of 
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the normal neutrosophic weighted Bonferroni mean (NNWBM) operator and the normal neutrosophic 
weighted geometric Bonferroni mean (NNWGBM) operator. In section 5, a numerical example is given 
to verify the proposed approach and to prove its effectiveness and practicality. In Section 6, we 
conclude the paper and give some remarks. 

2. Preliminaries 

2.1 The normal fuzzy set and normal intuitionistic fuzzy set 

Definition 1 [17], Let X be a real number set. A is denoted as A = (a, a) . If its membership function 
satisfies: 

A(x) = e a (ct > 0) (1) 

then A is called a normal fuzzy number. The set of all normal fuzzy numbers is denoted as N . 

Definition 2 [18,19]. Suppose X is an ordinary finite non-empty set and (a, a) eiV,A = ^(fl,cr),// A ,v A ) 

is a normal intuitionistic fuzzy number (NIFN) when its membership function is expressed as: 

_ ( 5 = 5) 2 

jU A (x) = n A e a ,xeX, (2) 

and its non-membership function is expressed as: 

_(AZ1) 2 

V',(.V) l-a-v',)*' ,xgX. (3) 

where 0 < x) < 1, 0 < v A ( x) < 1 and 0 < ju A + v A < 1 . When //, = I and v A =0, the NIFN 

will become a NFN. Compared to NFNs, the NIFN adds the non-membership function, which 
expresses the degree of not belonging to {a, a ) . Moreover, n A (x) = 1 - jU A (x) - v A (x) shows the degree 
of hesitance. The set of NIFNs is denoted by NIFNS. 

2.2 The neutrosophic set 

Definition 3 [8]. Let X be a universe of discourse, with a generic element in X denoted by x. A 
neutrosophic number A in X is expressed as: 

Mx) = (x I (T a (x), / A (x),F a (x))) (4) 

where, T A (x) is the truth-membership function, I A (x) is the indeterminacy-membership function, and 

F a (x) is the falsity-membership function. T A (x), I A (x) and F A (x) are real standard or nonstandard 

subsets of ](T,l + [. 

There is no restriction on the sum of T A (x) , I A (x) and F A (x) , so 0” < T A (x) + 1 A (x) + F A (x) < 3 + . 

Definition 4 [9]. Let X be a universe of discourse, with a generic element in X denoted by x. A single 
valued neutrosophic number A in X is 

A(x) = (x I ( T a (x),/ a (x), F a (x))) ( 5 ) 

where L, (x) is the truth-membership function, / A (x) is the indeterminacy-membership function, and 
F a (x) is the falsity-membership function. For each point x in X, we have T A (x), I A (x), F A (x) e [0,1] , 
and 0 <T a (x) + 1 A (x) + F A (x) < 3 . 
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2.3 The normal neutrosophic set 

Definition 5. Suppose X is a universe of discourse, with a generic element in X denoted by x, and 
(a, cr) e N , A normal neutrosophic number A in X is expressed as: 



AO) = (x I (ci,cr\(T A (x)J a (x),F a (x))) (6) 

Where the truth-membership function T A {x ) satisfies: 

T A (x) = T A e a ,xgX, 

the indeterminacy-membership function I A (x) satisfies: 

I A (x) = l-(l-I A )e * ,x e X. 
and the falsity-membership function F A (x) satisfies: 

-o=o 2 

F A (x) = l-(l-F A )e ° ,xgX. 

For each point x in X, we have T A (x), I A (x), F A (x) e [0, 1] , and 0 <T A (x) + 1 A (x) + F A (x)< 3 . The set 
of all normal neutrosophic numbers is denoted as R . 

Examplel. The service life of the lamp bulb obeys the normal distribution, the normal fuzzy number is 
IV (1000 ,302 ) . The experts evaluate whether the service life conforms to the normal distribution. At last, 
the experts give the evaluation values: the degree of result in range (1000 ,302 ) is 0.6; the degree of 
result not in range (1000 ,302) is 0.2; and the degree of hesitance is 0.2. So, the final evaluation result 

about the service life of the lamp bulb is A = ( (l 000,302), (0.6,0.2,0.2)) . 



Definition 6. Let =/(a 1 ,<7 1 ),(7J,/ 1 ,Fj)) and a 2 =((a 2 ,cr 2 ),(7^,/ 2 ,F 2 )) be two NNNs, then the 



Euclidean distance between a t and a, is defined as follows: 



d (x, y) = — 



J l{ 2 + T l L ~ ^ - F l L ) a i~{ 2 + T 2 ~ ~ F 2) a ij 



4 ' + ^ [ ( 2 + t; ^ - /,* - f; ^ ) <7, - ( 2 + T 2 l - 1 \ - F^ ) a 2 ] : 1 



(7) 



According to the operational laws defined by Wang et al. [22], we can give the following definition. 
Definition 7.Let 5j =^(o 1 ,(J 1 ),(7],/ 1 ,F 1 )) and a 2 = {^a 2 ,a 2 ),(T- L ,I 2 ,F 2 'fj be two NNNs, then the 
operational rules are defined as follows. 

(1) ®a 2 = +a 2 ,<T, +cr 2 ),(7j +T 2 -T i T 2 JJ 2 ,F l F 1 )) (8) 



2 2 

cr, cr. 



(2) d l ®a 1 =Ua l a 2 ,a l a 2 ]~Y + -f),(T l T 2 ,I i +I 2 -I l I 2 ,F i +F 2 -F l F 2 )) 

(3) A>0 

(4) a? = ((a 2 , A 1/2 <-' CT| ), (?;M - (1 - /, ) A ,1 - (1 - F, f )) A > 0 



(9) 



( 10 ) 



(11) 
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Theorem 1 . Let a l =([a l ,cr l ),(T 1 ,I 1 ,F l '^ and a 2 =((a 2 ,<7 2 ),(T,,/ 2 ,F 2 ))be two NNNs, and/7, 77,, 77, >0, 
then we have 



(DSj ®a n =tL ®a l 


(12) 


(2) ® = <7, ® a l 


(13) 


(3) 77(0, ©a 2 ) = 77 a x ®rja 2 


(14) 


(4) 77, a, ©77,^ =(77, + rj 2 )a t 


(15) 


(5) a q ® a 2 = (a, ® a 2 ) n 


(16) 


(6)3* Vo* = d 1 7l+ " 2 


(17) 



Definition 8. Let a k = ({a i ,a t 'j,(T k ,l k ,F t he a NNN, and then its score function is 

s 1 (d lc )=a k (2 + T k -I k -F k ), s 2 (a k ) = a k (2 + T k - I k - F k ) (18) 

and its accuracy function is 

h t Q k ) = a k (2 + T k -I k + F k ), h 2 {a k ) = ( J k (2 + T k -I k + F k ) (19) 

Definition 9. Leta] anda, be two NNNs, the values of score 

functions of^and d 2 are s, (a, ) , s 2 (a , ) and .v, (a 2 ) ,s 2 (a 2 ), and the values of accuracy functions offlj 
and cl, are,h x {d x ) , h 2 (a { ) and h x (d 2 ) ,h 2 (a 2 ) , respectively. Then, there will be: 

( 1) If jj(a[) > S[(a 2 ) 1 thenflj > a 2 ; 

(2) If i j («!) = .?! (a 2 ) , then 

® If h x (flj ) > /i[ (a 2 ) , then d x > a 2 ; 

© If /Z[ («j ) = h x (d 2 ) , then 

(i) If s 2 (d \ ) < s 2 (a 2 ) , then a x >a 2 \ 

(ii) If ,s 2 (fli) = s 2 (d 2 ) , then 

(a) If h 2 (a x ) < h 2 (a 2 ) , then d x > a 2 ; 

(b) If h 2 (a, ) = h 2 ( a 2 ) , then a l = a 2 . 

3. Normal neutrosophic Bonferroni mean operators 

3.1 NNBM and NNWBM operators 

Bonferroni [13] firstly introduced the Bonferroni mean (BM) which can provide the aggregation 
between the max and min operators and the logical “or” and “and” operators. However the Bonferroni 
mean (BM) operator [13] has mostly been used in the situation where the input arguments are the 
non-negative real numbers. In this section, we will study the BM operator under the environments of 
NNNs. Based on Definition of BM [13], we define the Bonferroni mean operator of NNNs as follows: 

Definition 10[13]. Suppose p,q>0 and {a 1 ,a 2 ,...,a m } is a set of NNNs. The Bonferroni mean operator 



of NNNs is defined as 



NNBM p ' q (a 2 ,a 2 ,...,d m ) 



1 



1(1 






J 



( 20 ) 
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Theorem 2. Leta t . ={^a k ,a k ),{T k ,l k ,F k ^j (k = 1,2 be a set of NNNs, then, the result aggregated 



from Definition 10 will be still a NNN, and even 



NNBM p ’ q (d l ,d 2 ,...,d m ) = { 



— — Y 



p Q 

a; a. 



r \ 






Z p-l o-l / 2 2 , 2 2 \ , 

a. a 1 pa j a i +qa.cr. I 2 



V <*2 J <*i 



J 



sjp + q ( m - 1) 



1- 



m 

Q(l-7 ] p T q ) 



i = i 
V >*i 



i 7 

(,„-!) 



, 1 - 



1- 



m / 

n l-C-t 

i.i = 1 v 



2 

V i*2 



n(m-l) 



( 21 ) 



1- 



1- 



<,2=1 V 



V 



1 3 



i 3 



Proof. By the operational rules of the NNNs, we have 

f i A 



32 


of = 


a-,P 


33 




V 


34 






35 




4 i 


36 


5 J = 


aj,? 2 


37 




V 


38 






39 


and 




40 




/, 


41 


=( 


42 




V 


43 






44 


then 




45 






4 6 


m 




47 




’®a« 


48 


U=l 

i*j 





,(7V\i-(i-/,.) p ,i-(i-.f;) p ) 



a? a % a?- 1 ay 1 (pa]af + qafcrj f j , ( T>Tf , 1 - (l • - ) P (l - / , )“ , 1 - (l - F i ) P (l - F . )" jj 



Z P <7 0-1 0 - 1 / 2 2 , 2 2\ 9 

<2,<2,’ZX a j \P a i CJ i +< l a i Cr j) 2 

<>2=1 <,2=1 

<*2 1*2 



<', 2=1 
V <*1 

and 



<, 2=1 

<*2 



<J=1 

<*2 



m( m — 1 ) 

i*2 






a, a 



ra(ra-l) /J= 

**7 



m 

Z p-l o-l/ 2 2 , 2 2\ 

a; aj \pa j (J i +qa i <j j \ 



i’j = 1 
**7 
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which completes the proof of the theorem 2. 

In the following, we will discuss some properties of NNBM operator as follows. 

Theorem 3. (Idempotency). 

Let{fl p a 2 ,...,d )n } be a set of NNNs, if all d k (k = 1,2,..., m) are equal, i.e.a^ = a (k = l,2,...,m) , 
then 



NNBM (d l ,d 2 ,...,a /ll ) = d 




which completes the proof of theorem 3. 

Theorem 4. (Commutativity). 

Let d[ (k = 1, 2, .... m) is any permutation of d k (k = 1, 2, .... m) .Then 
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NNBM ( Oj , a 2 , . . . , a m ) = NNBM (a[,a 2 ,...,a m ) 



Proof. Let NNBM (a, , d 2 , . . . , a m ) = 



NNBM p ’ q (a[,a' 2 ,...,a' m ) = 



1 






m ( m 


- 1 ) 


/ , a i 

ij= 1 


V 




i*j 






1 

p+q 


m 






2>7 






ij= 1 






i*j 


) 





Since {a[,a 2 ,.. a f n } is any permutation of {a l ,a 2 ,...,a m } , then we have fjif' . 



Thus, NNBM(d l ,d 2 ,...,d m ) = NNBM (d[,a 2 ,...,a m ) ■ 



which completes the proof of the theorem 4. 

Now we discuss some special cases of the NNBM by assigning different values to the 
parameters P> C 1 : 

(1) If r? =0,then 



NNBM p '° d m ) 



-S 

m , = i j 



( 22 ) 



which we call it the normal neutrosophic generalized mean (NNGM) operator. 

(2) If p = 1 and q = 0 , then 

| m 

NNBM^(a l ,a 2 ,...,a m ) = -Y J d, 

which we call it the normal neutrosophic mean (NNM) operator. 

(3) If P = 2 and q = 0 , then 

i 

f ] m \2 

NNBM 2,0 (5j , a 2 , . . . , a m ) = -15? 

V m i=i J 



(23) 



(24) 



which we call it the normal neutrosophic square mean (NNSM) operator. 
(4) If p = 1 and <7 = 1, then 



NNBM 1,1 (a 1 ,d 2 ,...,a m ) = 



Z at i 



i(m- 1) <;,j=i 



(25) 



which we call it the Normal neutrosophic interrelated square mean (NNISM) operator. 

The NNBM operator just considers the relationship of the aggregated arguments but ignores the 
importance of their weights. In the following, we will define another Bonferroni mean operator, the 
normal neutrosophic weighted Bonferroni mean (NNWBM) operator, to overcome the shortcoming. 



Definition 11. Let {d l ,d 2 ,...,a Itt } be a set of NNNs. The weighted Bonferroni mean operator of NNNs is 



defined as 




NNWBM p,q 



f J p Y 

( a 1 ,a 2 ,...,d m )= — — E {widi) P ®{wjdi) q 

m I m - 1 ) u = i 

V V ' i*J J 



Where w = (w 1 ,w 2 ,..., w m Y is the weight vector of NNNs. a k (k = 1,2,..., to) 0 <w k <\(k = l,2,...,w) 

m 

and Y, w k =1 - 

/c=l 

Theorem 5. Let a k =(^a k ,<J k ),{T k ,I k ,F k ^ (k = l,2---,w) be a set of the NNNs, then, the result 
aggregated based on the Definition! 1 will be still a NNN, and even 



7 7 2 (wiaiY ®{w j a J ) 

i[m — lj i,j=i 



7 7 X ( WidiY ) 

t(n -1] ij=i 



'L(w i aiY®{wja j y t (widiY l ®{wjaj)' p{wjCij ) (w,cr, )" + <?( W/fl,) (w 7 cry) 

*'J=1 iJ = 1 ' 

i*j / 

Jp + q (p+a Jln(n-l) 



n(i-(i-(i-r ( )r®(i-(i-r,)) wjS 

‘.j=i ' 



l- l- ,n i [i-(i-/r)'’(i-/?) < 'j a- i- n|i-(i-Fr) P (i-Fr) 4 j 

y v'*/ y 



The NNWBM operator has the following properties: 

Theorem 6. (Idempotency). 

Let {a,, a 2 ,...,a m } be a collection of NNNs, if all 5 k (k = 1,2,.. .,m)are equal, i.e. <5,. =5 



(A: = 1,2,..., to) , for all k, then 



NNWBM (a l ,a 2 ,...,a m ) = 5 



The proof of the Theorem 6 can be easily completed with the same way as the Theorem 3. 
Theorem 7. (Commutativity). 

Let a[ (k = 1, 2, ..., to) is any permutation of a k (k = 1, 2, ..., to) .Then 
NNWBM ( a 1 , a 2 , . . . , d m ) = NNWBM (a[,a 2 ,..„ a m ) 



The proof of the Theorem 7 can be easily completed with the same way as the Theorem 4. 




3.2 NNGBM and NNWGBM operators 

Definition 12. Suppose p,q > 0 and (a l ,a 2 ,...,a m ( be a set of NNNs. The geometric Bonferroni mean 
operator of the NNNs is defined as 
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NNGBM p ' q {a l ,a 1 ,...,a m ) = 



f 






t(m-l) 



p + q 



n {pdi + qaj) 



i,j = 1 



(28) 



Theorem 8. Let a k = ((a k ,cr k ),(T k ,I k ,F k ^j (k = l,2---,m) be a set of the NNNs, then, the result 
aggregated based on the Definition 12 will be still a NNN, and even 
NNGBM pq (a l ,a 2 ,...,a m ) = 



p + q ij=i 

i*j 



m . . 1 / 

n (pai + qay ,n(m ~ l \ 



p + q 



1 £ {PCTi + qCTj) 

i( m -l)uf(p ai + qaj ) 2 



ft ( Pdi + qa 



ij = 1 
i*j 



(29) 
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[i- n n (i ^ 

Proof. By the operational laws of the NNNs, we have 
pa i = ((pa, ,pa,),(l-(l-T,Y, //' , 

qZj = ((<?«; , q°j ) , (l - a - Tj ,I q ,,Ff)) 



and 



Pdi + qdj - ( ( Pai + qa p Pai + qa 



then 






fl (pdi + qdj)-{ 

ij = 1 
i*j 



fl ( Pai + qaj ). 

i,j=l 

i*j 



r 

X 

> 7 = 

V 1 */ 



(p<T, + <7<Ty) 



‘V; 1 (Pa/ + 9a;)' 



fl ( Pai + qaj ) 

*.7=1 

i*j 



fl (l-(l-7’,) P (l-7’i) )’l- n (l-/f/?),l- n ( 1 -FfFj) 

<, 4=1 ' ' i,j=i i,i = 1 






and 
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fl ( ptii+qtij ) 

i,j = 1 
\i*j 



n~ 1 ) 
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then 
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\i*j 
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which completes the proof of the theorem 8. 

The geometric Bonferroni mean operator of the NNNs has some properties as follows: 

Theorem 9. (Idempotency). 

Let {flj , a 2 , a m } be a set of the NNNs. If all a k (k = 1,2,..., m) are equal, i.e. a k = a (k = 1,2,..., m) , 
for all k, then 

AWG5M(aj,a 2 ,...,a m ) = a 

The proof of the Theorem 9 can be easily completed similar to the Theorem 3. 

Theorem 10. (Commutativity). 

Suppose (k = 1,2,..., m) is any permutation of a k (k = 1,2,..., m) .Then 



NNGBM (a v a 2 ,...,a m ) = NNGBM (a[, a 2 ,...,a m ) 

The proof of the Theorem 10 can be easily completed with the same way as the Theorem 4. 

Now we discuss some special cases of the NNGBM by assigning different values to the 
parameters P, q : 

(1) If <7 = 0, then 
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NNGBM p ’°(a i ,a 2 




( 30 ) 



which we call it the Normal neutrosophic generalized geometric mean (NNGGM) operator. 
(2) If p = 1 and q = 0 , then 

i 

NNGBM'' 0 {d 1 ,a 2 ,...,d m ) = [fl(ai) ) j m ( 31 ) 



which we call it the Normal neutrosophic geometric mean (NNGM) operator. 
(3) Up = 2 and 9 = 0, then 



NNGBM 20 (a 1 ,a 2 



1 

«,„) = |n(2a,)j 



(32) 



which we call it the Normal neutrosophic square geometric mean (NNSGM) operator. 
(4) If p = 1 and 9 = 1, then 



NNGBM ll (d 1 ,d 2 ,...,a m 



1 

'w(m-l) 



fl (di + dj) 

i,j = 1 



(33) 



\>*i J 

which we call it the Normal neutrosophic interrelated square geometric mean (NNISGM) 
operator. 

Similar to the NNBM operator, the NNGBM operator also just considers the interrelationship of 
the input arguments and ignores their own importance. In the following, we will extend the NNGBM to 
the normal neutrosophic weighted Bonferroni mean (NNWGBM) operator which can not only 
considers the interrelationship but also takes the weights into account. 



Definition 13. Let{d l ,d 2 ,...,a m } be a set of NNNs. The weighted geometric Bonferroni mean operator 



of the NNNs will be defined as: 



NNWGBM p ' q (d 1 ,d 2 ,..„ d m ) = 



p + q 






n \ p{d) w ' + q(dj)" 



<i(m- l) 



(34) 



Where w = (w l ,w 2 ,...,w m f is the weight vector of NNNs, a k (k = 1,2,..., m) ,0<w t < 1 (k = l,2,...,m) 

m 

and =1. 

k = 1 

Theorem 11. Let d k ={[a k ,a k ), (T k , 1 k ,F k ^j (k = l,2---,m) be a set of the NNNs, then, the result 
aggregated based on the Definition 13 will be still a NNN, and even 



NNWGBM (a, , fl 2 , ..., d rn ) = I 



— n( 

1 + q i,t=i V 
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The weighted geometric Bonferroni mean of the NNNs has some properties as follow. 
Theorem 12. (Idempotency). 



Let {aj ,a 2 ,...,a m } be a set of NNNs, if all a k ( k = 1,2,..., to) are equal, i.e. a k = a (k = l,2,...,m) , for 



all k, then 

NNWPG(a l ,a 2 ,...,a m ) = a 

The proof of the Theoreml2 can be easily completed with the same way as the Theorem 3. 
Theorem 13 (Commutativity). 

Let a[ ( k = 1, 2, ..., to) be any permutation of a k (k = 1,2,..., to) ■ Then 



NNGBM (a,, a 2 , ) = NNGBM (a[,a 2 ,...,a m ) 

The proof of the Theorem 13 can be easily completed similar to Theorem 4. 

4. A multiple attribute decision making method on the basis of NNWBM and 
NNWGBM operator 

In this section, we will apply the normal neutrosophic weighted geometric Bonferroni 
mean(NNWBM) operator (or NNWGBM) to solve the multiple attribute decision making problems on 
the basis of the NNNs. 

For a multiple attribute decision making problem, suppose A = {A l ,A 2 ,---, A„} is the set of the 

alternatives, and C = {C 1 ,C 2 ,---,C n } is the set of the attributes. Suppose each attributes are independent, 
and the evaluation value of the alternative A on the condition of the attribute C i is 
= ^(cl.ct), (77, Lj, Fy)^ , which is presented by the form of the NNN, where T . / . /■] e[0,l] and 

n 

T i} + Ijj + F tj < 3 . The weight vector of the attribute is w= (w 1 ,w 2 ,---,w n ) , which Wj e [0,1], /.w, =1. 

3=1 

Then, we use the normal neutrosophic weighted geometric Bonferroni mean(NNWBM) operator 
(or NNWGBM ) to develop a method to deal with the multiple attribute decision making problems as 
follows. 

Step 1. Normalize the decision matrix. 

Because there are two types of attribute, i.e., the benefit type and the cost type, we firstly convert 
the different types to the same one. So, the decision matrix of normal neutrosophic variables 

D = («, ; ) mx „ will be converted to the standardized matrix I) = (a (/ ) mx „ 




(36) 



For the benefit type: cl = 



max ; (a. ) max ( (<x ) 



For the cost type: 



\ a» max , ) a. 



Step 2. Calculate the comprehensive evaluation values of the alternatives based on the NNWBM 
operator (or NNWGBM). (generally, we can take p = q = l) 



it - NNWBM p,q (a fl , a i2 ,. . .,a in ) — / — — Y (w ja,ij\ ® f } 
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where i = 1,2, . 

Step 3. Calculate the score value of each comprehensive evaluation value by equation (18). 



Step 4. Rank all the al ternati ves { /t, , , ■ • • , /l m } and select the most desirable one(s) according to the 



definition 9. 

Step 5. End 

5. The numerical example 

In this section, based on NNWBM operator (NNWGBM), a numerical example is given to verify 
the proposed approach. 

There is a company which is planning to invest some money to an industry (cited from [10]). 
There are four alternative companies to be chosen, including: (1) A 1 is a car company; (2) A 9 is a food 
company; (3) A 3 is a computer company; (4) A 4 is an arms company. There are three evaluation 
attributes, including: (1) C t is the risk; (2) C 2 is the growth; (3) C 3 is the environment. We can know the 
attributes C) and C, are benefit criteria, and the type of C 3 is cost. The weight vector of the attributes 
is co = (0.35,0.25,0.4) . The final evaluation outcomes are expressed by the NNNs, and shown in table 1. 
Table 1 The evaluation values of four alternatives with respect to the three attributes 



Cl C2 C3 



Al 


<(3,0. 4), (0.4,0. 2,0. 3)> 


<(7, 0.6), (0.4,0. 1,0. 2)> 


<(5,0.4), (0.7,0. 2,0. 4)> 


A2 


<(4,0.2), (0. 6,0. 1,0.2)> 


<(8,0.4), (0.6, 0.1, 0.2)> 


<(6, 0.7), (0.3,0. 5,0. 8)> 


A3 


<(3.5, 0.3), (0.3, 0.2,0. 3)> 


<(6,0.2),(0.5,0.2,0.3)> 


<(5.5,0.6),(0.4,0.2,0.7)> 


A4 


<(5,0.5), (0.7,0. 1 ,0.2)> 


<(7,0.5),(0.6,0.1,0.1)> 


<(4.5,0.5),(0.6,0.3,0.8)> 



5.1 Procedure of decision making method based on the NNWBM operator. 

(1) Normalize the decision matrix 

Since Cl and C2 are benefit attributes, and C3 is a cost attribute, we utilize the formulas (36) and 
(37) to obtain the standardized decision matrix, which is shown in Table 2. 



Table 2. The standardized decision matrix 





Cl 


C2 


C3 


Al 


<(0.6,0. 1067),(0.4,0.2,0.3)> 


<(0.875,0. 0875), (0. 4,0. 1,0.2)> 


<(0.9,0. 0475), (0.4,0. 8, 0.7)> 


A2 


<(0. 8,0. 02), (0.6,0. 1,0.2)> 


<(1,0. 0333), (0. 6,0. 1,0.2)> 


<(0.75,0. 1167),(0.8,0.5,0.3)> 


A3 


<(0.7,0.05 14), (0. 3,0.2, 0.3)> 


<(0.75,0.01 11), (0.5,0.2,0.3)> 


<(0.8 18,0. 0935), (0. 7,0. 8,0.4)> 


A4 


<( 1 ,0. 1 ),(0.7,0. 1 ,0.2)> 


<(0.875,0. 0595), (0. 6,0. 1,0.1)> 


<( 1,0. 0794), (0.8,0. 7,0. 6)> 



(2) Calculate the comprehensive evaluation value of each alternative by formula (38). (suppose p=q=l). 



5, = ((0.1827, 0.0208), (0.5704, 0.7848, 0.8133)) 
a 2 = ((0.1954,0.0169), (0.6111,0.7084,0.7258)) 

03 = ((0.1761, 0.0143), (0.6232,0.8102,0.7834)) 
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a A = ((0.225 1,0.0190), (0.5770, 0.7419, 0.7535)) 

(3) Calculate the score function by formula (18). 

= 0.1776, sj(a 2 ) = 0.2299,5,(^3) =0.18 13, s, ((34) = 0.2435 

(4) Rank all of the alternatives and choose the most desirable one by the score function. 

According to the score function 5,(0 ) , the ranking is A 4 >- A, >- A, y A l . 

Thus, the best alternative is A 4 . 

5.2 Procedure of decision making method based on the NNWGBM operator. 

(l)Normalize the decision matrix 

Since Cl and C2 are benefit attributes, and C3 is a cost criterion, we use the formulas (36) and (37) 
to get the standardized decision matrix, which can be shown in Table 3. 



Table 3. The standardized decision matrix 





Cl 


C2 


C3 


A1 


<0.6,0. 1067), (0.4,0.2,0.3)> 


<0.875,0. 0875), (0. 4,0. 1,0.2)> 


<0.9,0. 0475), (0.4,0. 8, 0.7)> 


A2 


<(0. 8,0. 02), (0.6,0. 1,0.2)> 


< 1,0. 0333), (0. 6,0. 1,0.2)> 


<(0.75,0. 1167),(0.8,0.5,0.3)> 


A3 


<0.7,0.05 14), (0. 3,0.2, 0.3)> 


<0.75,0.01 11), (0.5,0.2,0.3)> 


<0.8 18,0. 0935), (0. 7,0. 8,0.4)> 


A4 


<1,0.1), (0.7, 0.1, 0.2)> 


<0.875,0. 0595), (0. 6,0. 1,0.1)> 


<( 1,0. 0794), (0.8,0. 7,0. 6)> 



(2) Calculate the comprehensive evaluation value of each alternative by formula (39). (suppose p=q=l). 



a t = ((0.6783, 0.0302), (0.8143, 0.0917, 0.1101)) 
a 2 = ((0.6850, 0.0224), (0.8556, 0.0517, 0.0596)) 
ch = ((0.6748, 0.0207), (0.8567,0.1050,0.0892)) 

a 4 = ((0.7032, 0.0240), (0.8372, 0.0643, 0.0744)) 

(3) Calculate the score function by formula (18). 

si(5i ) = 1-7721, 5i (o 2 ) = 1-8798, 5, (a 3 ) =1-7968, 5, (d 4 ) =1-8977 

(4) Rank all of the alternatives and choose the most desirable one by the score function. 

According to the score function 5, (a .) , the ranking is A, >- A, >- A, A A, . 

Thus , the best alternative is A, . 

5.3 Analysis the effect of the factor 7 

In order to demonstrate the influence of the parameter 7 on decision making results of this 
example, we use the different values/ 7 ’ 7 in NNWBM or NNWGBM operator in step 4 to rank the 
alternatives. The ranking results are shown in Table 4 and Table 5. 

Table 4 Ordering of the alternatives by utilizing the different Pi 7 in NNWBM operator 
P>7 



Score values 5 1 (« / ) 



Ranking 





p=Q,q = l 



5,(5,) = 0.1292,5,(5 2 )= 0.1457, 
5,(a 3 ) = 0.1235, 5,(5 4 ) = 0.1648 



A 4 >- A, >- A i > A 3 



p = 0,q=2 



5,(5,) =0.2357, 5,(5 2 ) =0.2517, 
5,(5 3 ) = 0.2206, 5,(5 4 ) = 0.3003 



A 4 >- A, y A, >- A 3 



p=0,q = l0 



5, (5, ) = 0.4494, 5, (5, ) = 0.4448, 
5,(a 3 ) = 0.4163, 5,(5 4 ) = 0.5902 



A 4 >- A, >- A, X A, 



p = !,<? = 0 



5,(5,) = 0.0879, ,v, (a 2 ) = 0.1338, 
5,(5 3 ) = 0.1002, .s', (5 4 ) = 0.1486 



A 4 >- A, >- A 3 >- A, 



p = 2,q = 0 



.s, (a, ) = 0.1569,5,(5,) = 0.2394, 
5,(5 3 ) = 0.1756,5,(5 4 ) = 0.2679 



A 4 >- A, >- A 3 >- A, 



p = l0,q = 0 



5, (5,) = 0.2772, 5, (5,) =0.4 181, 
5,(5 3 ) = 0.3 158, 5, (5 4 ) = 0.5 1 66 



A 4 >-A 2 >A 3 > A, 



p = 2,q = l 



5, (5, ) = 0.2087, 5, (5, ) = 0.2824, 
5,(5 3 ) = 0.2179, 5,(5 4 ) = 0.3061 



A 4 >- A, >- A 3 >- A, 



p = 10,g = l 



5,(5,) =0.2853, 5,(5,) =0.4185 
5,(5 3 ) = 0.3181, 5,(5 4 ) =0.5075 



A 4 >-A 2 >A 3 >- A, 



p=l,q=l 



5,(5,) = 0.1776,5,(5,) =0.2299 
5, (5 3 ) = 0. 1 8 1 3, 5, (5 4 ) = 0.24 35 



p =\,q = 2 



5,(5,) = 0.2408,5,(5,) = 0.2848 
5,(5 3 ) = 0.2350, 5,(5 4 ) = 0.3161 



A 4 > A, > A x ^ A^ 



p =l,g =10 



5,(5,) = 0.4270, 5,(5,) = 0.4328, 
5,(5 3 ) = 0.3942, 5,(5 4 ) = 0.5552 



A 4 > A,_ > A { > A 3 



Table 5 Ordering of the alternatives by utilizing the different P' c ! in NNWGBM operator 
P’ c l Score values 5,(5, ) Ranking 



p = 0,q = l 



5,(5,) = 2.5272, 5,(5,) = 2.6809, 
5,(5 3 ) = 2.5 721, 5, (5 4 ) = 2.65 87 



A, >- A 4 x A 3 x A, 



p = 0,q = 2 



5,(5,) = 1.6499,5,(5,) =1.8015, 
5,(5 3 ) = 1.7067 5, (5 4 ) = 1.72 47 



A, >- A 4 >- A 3 >- A, 
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p =0,q = 10 



.s', (5, ) = 0.6120, .s, (a 2 ) = 0.7106, 
5,(5 3 ) = 0.6521, .5,(54) = 0.6168 



A, >- A, >- A 4 >- A, 



p = l,q = 0 



5,(5,) = 2.5272, 5,(5,) = 2.6809, 
5,(5 3 ) = 2.5 721, 5, (5 4 ) = 2.65 8 7 



\ >- A 4 >- A 3 >- A, 



p = 2,q = 0 



5,(5,) = 1.6499,5,(5,) =1.8015, 
5,(5 3 ) = 1.7067, 5,(5 4 ) = 1.7247 



A, ^ A 4 x A 3 x A, 



p = 10,g = 0 



5,(5,) = 0.6121,5,(5,) = 0.7106, 
5,(5 3 ) = 0.6521, 5, (5 4 ) = 0.6168 



A 1 >~ A 3 >~ A 4 >~ A, 



P = 2, q = l 



5,(5,) = 1.3853, 5, (5 2 ) = 1.5031, 
5,(5 3 ) = 1.4232, 5,(5 4 ) = 1.4891 



A, >- A 4 >- A 3 >- A, 



p = l0,q = l 



5,(5,) = 0.6081,5,(5,) = 0.7005, 
5,(5 3 ) = 0. 6433,5, (5 4 ) = 0.6295 



A 1 >- A 3 >~ A 4 >- A, 



p = l,q = l 



5,(5,) = 1.7721, 5, (5,) = 1.8798 
5,(5 3 ) = 1.7968, 5, (5 4 ) = 1.8977 






p = l,q = 2 



5,(5,) = 1.4110, 5,(5 2 ) =1.4988 
5, (5 3 ) = 1 .4304, 5, (5 4 ) = 1 .5024 



A, >- A, y A, >- A, 



p =l,q =10 



5,(5,) = 0.6236, 5,(5,) = 0.6972, 
5,(5 3 ) = 0.6486, 5,(5 4 ) = 0.6370 



A, >- A 3 >- A 4 x A, 



As we can see from Table 4, the ordering of the alternatives may be different for the different 
values of P^l in NNWBA operator. But the best alternative is the same one A 4 . In the Table 5, the 

ordering of the alternatives also may be different for the different values of 74 7 . The best alternative is 
A,orA 4 .In practical applications, we generally adopt the values of the two parameters as/? = <7 = 1, 
which are not only easy and intuitive but also fully capture the correlations between criteria. 

6. Conclusions 

The multiple attribute decision making method on the basis of normal neutrosophic variables has a 
wider application in many domains. The normal neutrosophic set (NNS) will be more appropriate to 
deal with the incompleteness, indeterminacy and inconsistency of the decision-making information, 
and the Bonferroni mean (BM) operator can consider the interrelationships between the input 
arguments. So, in this paper, we proposed two aggregation operators called the normal neutrosophic 
Bonferroni mean (NNBM) operator and the normal neutrosophic geometric Bonferroni mean 
(NNGBM) operator for aggregating the information expressed by the normal neutrosophic numbers. 
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We studied some properties of them and discuss their some special cases. For the situations in which 
the input arguments have different weight, we then developed the normal neutrosophic weighted 
Bonferroni mean (NNWBM) operator and the normal neutrosophic weighted geometric Bonferroni 
mean (NNWGBM) operator, on the basis of which we propose two procedures for multiple attribute 
decision making under the environments where the information is expressed by the NNNs. Moreover, 
we use the NNWBM operator and NNWGBM operator to aggregate the evaluation information of 
alternatives, so the decision makers can get the desirable alternative according with their interest and 
the practical need by changing the values of P’Q , which makes the results of the proposed multiple 
attribute decision making method more flexible and reliable. In the further research, the study about the 
applications of the new decision making method is necessary and significative because the applications 
of the normal distribution is widely distributed in many domain in the uncertain environment. 
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